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Verfahren und Anordnung zur Detektion und -zur Messung 
der Phase von periodischen Biosignalen 

5 Die Erfindung betrifft ein Verfahren und eine Anordnung zur 
echtzeitfahigen zuverlassigen Detektion und Messung der Phase 
periodischer physiologischer GroSen oder Biosignale. 

Im Stand der Technik sind Verfahren bekannt, die zur 
10 Bestimmung der Phase ein uber den zeitlichen Verlauf des 
Biosignals gleitendes oder sequentiell angeordnetes 
Analysef enster verwenden. Auf den im Analysef enster liegenden 
Signalabschnitt werden auf der Fourier-Trans formation 
basierende Methoden oder deskriptive Statistiken angewandt. 
15 So werden beispielsweise in der Perimetrie periodisch 
auf leuchtende Lichtmarken definierter Intensitat und 
hinreichend hoher Frequenz (uber etwa 4Hz) beriutzt, um die 
Funktionsfahigkeit des visuellen Systems zu uberprufen. Fur 
den Funktionstest wird- das Elektroenzephalogramm (EEG) erfasst 
20 "und die Reizantwort auf den visuellen Stimulus bezuglich der 
Amplituden und der Phase analysiert. Die Phase der Reizantwort 
ist ieiner der entscheidenden diagnostischen Parameter in der 
Funktionsdiagnostik . 

2 5 Bei den bisherigen Verfahren ist nachteilig, dass die 

statistische Unsicherheit der Detektion bzw. die Ungenauigkeit . 
der Messung sehr hoch ist. Die Unsicherheit und die 
Ungenauigkeit ergeben sich aus der Signaltheorie als Folge von 
und im Zusammenhang mit der Lange des Analysef ensters . Die 

3 0 Theorie besagt, dass mit abnehmender Lange des Analysef ensters 

die statistische Unsicherheit und damit die Ungenauigkeit 
zunehmen, was in der praktischen Signalanalyse auch 
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hinreichend bewiesen und bekannt ist. Far ein statistisch 
besseres Ergebnis miisste zunachst die Fensterlange vergroSert 
werden. Aus der Physiologie jedoch ist bekannt, dass sich die 
Phase relativ schnell andern kann und. diese Anderungen auch 
diagnostisch relevant sind. Bei einem langen Analysefenster. 
geht die wertvolle Information uber die Phasenanderung 
verloren und die statist! sche Unsicherheit des Messergebnisses 
nimmt infolge der Anderungen nicht zwangslaufig ab. 

Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren und 
eine Anordnung anzugeben, . mit denen es moglich ist, einen 
kausalen Phasengang in periodischen Biosignalen mit einer ■ 
gegeniiber herkommlichen Verfahren besseren Zuverlassigkeit und 
h6heren Geschwindigkeit bei gleichzeitig reduziertem 
Rechenaufwand zu detektieren und zu messen. 

Erfindungsgemafc gelingt die L6sung der Aufgabe dadurch, dass 
periodische Biosignale entsprechend ihrem physikalischen und 
physiologischen Ursprung erfasst werden, dass parallel zum 
analysierten biologischen System ein - Zustandsbeobachter 
aufgestellt wird, dass die. Ausgangsgr6Sen des biologischen 
Systems und des Beobachters mit Hilfe eines Kalman-Fi Iters 
ausgewertet werden und zur Bestimmung der Phase genutzt 
werden . 

Bei dem erf indungsgemSfien Verfahren wird die Phase eines 
periodischen Biosignals ermittelt. und fur 

funktionsdiagnostische Zwecke genutzt. So ist beispielsweise 
eirie verlangerte Phase im Vergleich mit der gesunder 
Versuchsobjekte ein wichtiger Hinweis auf funktionelle 
Probleme des Untersuchten biologischen Systems. 
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Bei der erf indungsgem&Sen * Anordnung ist parallel zum 
untersuchten biologischen System, das mit einem Zustandsmodell 
nachgebildet wird, ein Zustandsbeobachter angeordnet, der 
entsprechend dem Systemmodell die Zustandsgrofie Phase 
basierend auf dem Kalman-Filter schatzt. 

Vom Vorteil ist hierbei, dass die Schatzung der Phase 
kontinuierlich stattfinden kann und kein gleitendes .oder 
sequentielL angewandtes Analysef enster notwendig ist. Damit 
wird die Analyse der zeitlichen Phasenanderungen erst moglich. 
Im Gegensatz zum irelativ komplizierten theoretischen 
Hintergrund dieses Phasenschatzers ist die praktische 
Umsetzung einfach. Im Vergleich zu herkommlichen Verfahren 
benotigt sie wesentlich weniger Rechenlei- stung, so dass eine 
Phasenschatzung in Echtzeit moglich ist. 

Die Erfindung wird im folgenden anhand der theoretischen 
Herleitung und eines Ausfuhrungsbeispieles naher erlautert. In 
den zugehorigen Zeichnungen zeigen: 

Figur 1 - Blockdiagramm des Beobachterkonzeptes 
Figur 2 - Systemmodell 

Figur 3 - Darstellung des Zustandsbeobachters zur Messung 
der Phase in periodischen Biosignalen 

Figur 4 - Verlauf der geschatzteh Phase fur eine 

Harmonische der Frequenz 8Hz und Phase 2rad fur 



Figut 5 



die statischen Kalman-Faktoren 2 und 20 
Verlauf der geschatzter^ Phase far eine 
verrauschte Harmonische der Frequenz 8HZ und 
Phase 2rad mit einem SNR von OdB (unten) und 
dynamischen Kalman-Faktor (oben) 
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Figur 6 - Phasenschatzung des Signals wie in Fig. 5 mit ' 

statischen Kalman-F^ktoren 
.Figur 7 - Ergebnisse der Phasenschatzung (rechts) 
an real en Signalen (links) 

Figur 8 - additive Uberlagerung einer Harmonischen einer 

Frequenz von 8Hz mit dem Rauschen und gezielter 
Verstimmung der ■ Analysef requenz (oben) und der 
Phasenverlauf mit Anstieg (linten) 

Das biologische System, das ein periodisches Biosignal 
produziert bzw. auf ein periodisches Eingangs signal antwortet, 
ist in Figur 1 als Zustandsmodell „reales System" dargestellt . 
Die Zustandsgleichungen (1) und (2) beschreiben dieses System 
(Fettgedruckte groSe Buchstaben stehen far Matrizen, kleine 
fur Vektoren) : 



x(t) = A-x(t) + B-u(t); x(0) = x(t 0 ) 
y(t) = C'X(t). 



(1) 
(2) 



Fur weitere Betrachtungen nehmen wir ein additives 
Signalmodell an, welches eine harmonische Schwingung und 
weifies normalverteiltes Rauschen summiert: 



y(t) = y.sin(cot + (p(t)) + r D (t) 



(3) 



Das Ziel • besteht darin, ein Systemmodell zu konstruieren, 
dessen Variable x(t) die Phase <p(t) des zu untersuchenden 
Signals y,(t) reprasentiert . Die Phase kann nicht direkt 
gemessen werden, da sie Argument einer trigonometrischen 
Funktion ist. Daher wird eine Hilf skonstruktion benotigt. Eine 
solche Konstruktion ist ein Zustandsbeobachter , der parallel 
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zum « untersuchten System, angeordnet ist. Der Beobachter schatzt 
die . Zustands variable mit Hilfe der Minimierung einer 
. Fehler funktion, welche . die Ausgange des realen Systems und des 
Beobachters vergleicht . Auf diese Weise kann nach 
5 abgeschldssener Fehlerminimierung die' Zustandsvariable Phase 
direkt gemessen werden. 



Das Blockdiagramm des Beobachterkonzeptes ist in Figur 1 
dargestellt. Da x(t) nicht direkt gemessen werdeii kann, wird 
10 Xra(t) im Beobachter geschatzt. • Die innere Schleife im 
Beobachter minimi ert den Fehler von y m (t) bezuglich y(t) mit 
• Hilfe der Korrekturmatrix K. Ftir den Beobachter ergeben sich 
dann die Zustandsgleichungen (4) und (5) . 

15 x M (t) = Ax M (t)+Bu(t) + K[y(t)-y M (t)} / (4) 

y M (t) = C-x^(t). (5) 

Aus (4) und (5) folgt: 

2 0 x M (t) = (A - KC) - x M (t) + B • u(t) + K ■ y(t) . ( 6 ) 

Eg wird davon ausgegangen, dass beide Systeme unterschiedliche 
Anfangsbedingungen haben. Daraus ■ ergibt sich der 
Beobachtungs fehler : 



25 



e(t) = x(t)-x M (t). (7 ) 

Der Beobachtungs fehler verschwindet iterativ mit Hilfe der 
Korrekturmatrix K, so dass 



30 
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e(t) =s 0 for t -> oo. . ( 8 ) 

Die Dynamik und die Stabilitat der Schatz.ung k6nnen mit der 

Differentialgleichung des Beobachtungs f ehlers (9) beschrieben 
werden : 

e(t) = x(t)-x M (t). (9) 

Durch Umstellung und weitere Zwischenschritte erhalten wir: 
e(t) = (A-KC).e(t). (10) 



Entsprechend dem Signalmodell (3) ist damit zu rechnen, dass 
das xintersuchte Signal durch Rauschen gestQrt ist. Urn den 
Einfluss des Rauschens zu reduzieren, wird ein Kalman-Filter 
eingesetzt. Das Rauschen berticksichtigend wird das System 
durch folgende Zustandsgleichungen beschrieben: 

Systemzustand: x(t) = A-x(t)+B-u(t) + M-r 8 (t) 
Systemausgang: y(t) = C-x(t) + r p (t) 

Beobachter : x M (t) = A • x M (t) + B • u(t) + K(t) • \y (t) - C x M (t)] , 
wobei 

K(t) ist die Korrekturmatrix, welche zu erreichen hat, dass 

e(t)=x(t)-x M (t)^0, 

e(t) ist der Beobachtungs fehler, 

r s (t)ist das Systemrauschen, und 

r p (t)ist das Prozessrauschen. ■ 



(11) 
(12) 
(13) 
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Urn die Herleituhg zu vereinf achen, wird angenommen, dass die 
Rauschanteile bre'itbandige GauSsche Nullmittelprozesse mit 
bekannten Kovarianzen sind: 

cov rp (t p t 2 ) = E{r p (t,) • r p T (t 2 )}= R p ( tl ) • 8(t, - t 2 ) 

cov r> (t p t 2 ) = EfcCt,) • r7(t 2 )}= R 8 ( tl ) • 8(t, - 1 2 ) ( 14 ) 

Die Rauschanteile sind voneinander unabhangig, also 

cov r p Jt 1 .t 2 ) = 0. (15) 

Far eine konsistente Schatzung . von x(t) muss die 
Fehlerleistung mit Hilfe der Matrix K(t) minimiert werden: 

E{e T (t) • e(t)}=E^ 2 (t) + e*(t) + • • • + e* (t)}= f (K(t))=Min. , (16) 

Unter Annabme der stochastischen Beziehungen bezuglich der 
Kovarianzen wird entsprechend dem Kalman-Filter eine geeignete 
Korrekturmatrix K(t) abgeleitet: 

K(t) = cov e (t)-C T -R; 1 (t). (17) 

Die Formel ftir die Fehlerkovarianz cov e (t) kann vom Kalman- 
Filter abgeleitet werden: 

cov e (t) = A.cov e (t)+^ (18) 
Schatzung der Phase: 

Das untersuchte Signal wird entsprechend (3) aus der Suinme 
einer Harmonischen und Rauschen modelliert: 
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y('0 = y nI (t)+i^(t) = y-sin(Q)t + (p(t))+i i( (t),Q)=27rf. (19) 

Die Phase ergibt sich aus der Dif ferentialgleichung: 

<ftt)=-a<p(t)+i;(t) a>0. (20) 

Daraus ergibt sich das Systemmodel in Figur 2. Die Phase kann 
nicht direkt gemessen werden. Daher liegt parallel zum System 
ein Beobachter, in dem ein -direkter Zugrif f auf die geschatzte 
Phase <pM(t) moglich ist. Allerdings der nichtlineare Anteil 
Yni(t) in (19) ist far das Beobachterkonzept ungunstig. Eine 
geeignete Linearisierung yi «p(t) , t) wird benStigt, um eine 
lineare Beziehung zwischen der Zustandsvariablen q>n(t) und dem 
Ausgang y M (t) entsprechend (5) herzustellen . Basierend auf der 
Taylor-Linearisierung kann der Beobachter . f olgendermaSen 
formuliert werden: 

<j) M (t) = -a • <p M (t) + K(t) • (y(t) - y M (t)) (21) 

y M (t) = yi(q> M (t).t), (22) 

Entsprechend (21) und (22) wird der Beobachter modelliert, wie 
in Figur 2 dargestellt. Im Ergebnis der Linearisierung im 
Arbeitspunkt qb, wird (22) mit (5) verknupf t . Daraus ergibt 
sich der Faktor C, welcher in (17) verwendet wird, um den 
Korrekturfaktor K(t) zu bestimmen: 

K(t) = cov e (t) • y • cos(Qt + <p B (t)) • R; l (t). (23) 
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Als Arbeitspunkt wird'die zu bestimmende Phase gewahlt 

<Pb(0=9m(0» (24) 

und die resultierende Dif ferentialgleichung fur die Phase ist 

9m (t) = -a • <p M (t) + cov e (t) • y • cos((»t + <p M (t)) • [y(t) - y . sin(cot + q> M (t))] • R" 1 (t). (25) 

Entsprechend (25) kann der Phasenschatzer modelliert werden, 
wie in Figur 3 dargestellt. 

Zur Phasenschatzung in y(t) muss die Fehlerkovarianz berechnet 
werden. Aus (18) folgt: 

cov e (t) = -2a • cov e (t) - y 2 • cos 2 (cot + 9m (t)) • cov 2 (t) • R; 1 (t) + R s (t). (26) 

Gleichung (26) ergibt eine einfache Losung,. falls 
hQherfrequente Anteile in der Fehlerkovarianz nicht 
beriicksichtigt werden. Basierend auf (27) 

cos 2 (cot + <p M (t)) = 1/2 + 008(2(01 + 2(^(1)) (2 7) 
kann (26) vereinfacht werden zu: 

cov e (t) = -2a.cov e (t)-l/2y 2 .R- I (t).cov 2 (t) + R s (t). (28) 

Bei geeigneter wahl des Parameters a in (28) . werden 
hochfrequente Anteile in Folge temporaler Integration 
unterdrtickt, also es liegt ein Tiefpassverhalten vor. Unter 
Berticksichtigung des Tiefpasses kann (25) vereinfacht werden: 

25-02 / 27.11.02 
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<p M (t) = -a-9 M ( t ) + cov e (t)y C Qs((0t + (p M (t))y(t)R^(t), (29) 

Damit wird der Beobachter, dargestellt in Figur 3, 
vereinfacht. Das in (29) vorgeschlagene System kann 
insbesondere zur Phasenschatzung von Harmonischen im Rauschen 
eingesetzt werden. 

In Figur 4 ist der Verlauf der geschatzten Phase fur eine 
Harmonische der Frequenz 8Hz und Phase 2rad far 
unterschiedliche Kalman-Faktoren dargestellt. Die Kalman- 
Faktoren sind statisch und betragen 2 bzw. 20. Wie der Grafik 
entnommen werden kann, mit ein'em hoheren Kalman-Faktor wird 
die Schatzung langsamer. Statische Kalman-Faktoren massen dort 
eingesetzt werden, wo der Zeitpunkt der Phasenanderung nicht 
bekannt ist. 

In Figur 5 ist der Verlauf (unterer Teil) der geschatzten 
Phase far eine verrauschte Harmonische der Frequenz 8Hz und 
Phase 2rad mit einem SNR ( Signal -Rausch- Verbal tnis) von OdB 
und dynamischen Kalman-Faktor (oberer Teil) dargestellt. Ist 
der Zeitpunkt der Phasenanderung bekannt, so kann der Kalman- 
Faktor so konstruiert werden, dass zunachst der Anderung 
gefolgt und anschliefiend die Varianz der Schatzung minimiert 
wird. Zum Vergleich ist in Figur 6 die Phasenschatzung des 
selben Signals wie in Figur 5 mit statischen Kalman-Faktoren 
dargestellt. 

In Figur 7 sind die Ergebnisse der Phasenschatzung (rechte 
Spalte der Grafik) an realen Signalen dargestellt. Zur 
Stimulation des visuellen Systems wurden Lichtpulsfolgen mit 
der Wiederholrate von 8 Pulsen pro Sekunde verwendet, wobei 
abwechselnd einer Ruhepause (Zeit von 0 bis 2 der Zeityeriaufe 
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der Signale oben links und • un ten links) die Lichtstimulation 
(Zeit 2 bis 5) folgte. Dargestellt ist der Verlauf des EEG 
(Elektroen'zephalogramm,' linke Spalte der Graf ik) von zwei 

okzipitalen Positionen nach 16-facher reizbezogener Mittelung. 

Beide Zeitverlaufe zeigen einen deutlichen Sprung der Phase 

nach Beginn der Lichtstimulation. 

Die Phasenschatzung wird problematisch bei stark verrauschten 
Signalen. Generell gilt zwar, dass die Phase robuster gegen 
Storungen ist als die Amplituden, wie dies schlieSlich auch in 
der Informationstechnik bekannt ist. Allerdings ist in diesem 
Grenzbereich zuhachst die Frage nach dem Vorhandensein - also 
der Detektion - einer kausalen Phase zu klaren, erst dann ware 
die Phase zu schatzen. In Figur 8 ist eine Harmonische der 
Frequenz 8Hz additiv dem Rauschen tiberlagert beginnend bei 
t=4s, wobei das SNR -10dB betrSgt (oberer Verlauf in der 
Grafik) . . im Zeitbereich ist die Harmonische beziiglich ihrer 
Amplituden nicht nachweisbar. Verwendet man den Phasenschatzer 
mit einer gezielten Verstimmung, hier mit einer Frequenz von 
7.8Hz, also urn 0.2Hz weniger als die Frequenz der 
Harmonischen, so ergibt sich im Falle einer kausalen Phase ein 
Anstieg von 1.2rad/s (unterer Verlauf der Grafik). Dieser 
Anstieg kann in Kombination mit einem Diskriminator direkt zur 
Detektion des Signals genutzt werden. 
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Bezugszeichenliste 



a - Systeinparameter, im Beobachtennodell wahlbar 



A, B, 




C, K 


Matrizen im Zustandsmodell eines Systems 


cp(t) 


Phase im Systemmodeli, zu schatzende GroSe 


<P M (t) 


Phase im Beobachtermodell , messbare Gr6£e 


r p (t) 


Pr o z e s s r aus chen 


r s (t) 


Systemrauscheh 


u(t) 


Eingangsvariable eines Systems im Zustandsmodell 


x(t) 


Zustandsvariable eines Systems im Zustandsmodell 


3£ M (t) 


Zustandsvariable. des Beobachters im Zustandsmodell 


y(t) 


Ausgangsvariable eines Systems im Zustandsmodell 


y M (t) 


Ausgangsvariable des Beobachters im Zustandsmodell 


yi(cp(-t) ,t) - 


Linearisierungsoperator fur Phase 



25 



30 
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Patentanspriiche 

Verfahren zur Detektion und. zur Messung der Phase von 
' pqriodischen Biosignalen, dadurch gekennzeichnet, dass 

a) das zu analysierende Signal, dessen Phase bestimmt 
werden soil, einem .Multiplizierer zugefuhrt wird, 
dessen Faktor zeitlich konstant oder veranderlich 
sein kann, 

b) das erhaltene Produkt einem weiteren Multiplizierer 
zugefuhrt wird, dessen Faktor eine trigdnometrische 
Funktion ist, deren Argument sich aus dem Produkt der 
untersuchten Frequenz mit der Zeit, addiert 'mit der 
gemessenen Phase, ergibt, wobei die Frequenz der 
trigonornetrischen Analysefunktion der Frequenz 
entspricht, bei der die Phase ermittelt werden soil, 
oder um einen bekannten Betrag von dieser Frequenz 
abweicht, 

c) die gemessene Phase einem dritten Multiplizierer 
zugeftthrt wird, in dem sie mit ' einem Faktor 
multipliziert wird, 

d) das unter b) erhaltene Produkt auf den positiven 
Eingang und das . unter c) erhaltene Produkt auf den 
negativen Eingang eines Dif f erenzbildners zugefuhrt 
werden, 

e) die. unter d) erhaltene Differenz einem zeitlichen 
Integrator zugefuhrt wird, an dessen Ausgang die zu 
ermitteinde gemessene Phase anliegt. 
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Anordnung zur Detektion und zur Messung der. Phase von 
periodischen Biosignalen, dadurch gekennzeichnet, dass 

a) parallel zum untersuehten biologischen System ein 
Zustandsbeobachter angeordnet wird, 

b) der Zustandsbeobachter die Funktionen nach Anspruch 1 
realisiert . 
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Zusammenf assung 

5 

1. Verfahren und Anordnung zur Detektion und zur Messung der 
Phase von periodischen Biosignalen 

2.1. Es ist ein Verfahren und eine Anordnung anzugeben, mit 
10 denen es moglich ist, einen kausalen Phasengang in 

periodischen Biosignalen ' mit einer gegenuber 

herkommlichen Verfahren besseren Zuverlassigkeit und 
hoheren Geschwindigkeit bei gleichzeitig reduziertem 
Rechenaufwand zu detektieren und zu messen. 

2.2. Parallel zum analysierten biologischen System wird ein 
Zustandsbeobachter *auf gestellt . Die AusgangsgroS.en des 
biologischen Systems und des Beobachters werden mit Hilfe 
eines Kalman-Filters ausgewertet und zur Bestimmung der 

20 Phase genutzt werden. 

2.3. Die ermittelte Phase eines periodischen Biosignals wird 
fur funktionsdiagnostische Zwecke genutzt. 

25 3 . Figur 1 
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